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Abstract 

We study the combined effect of nuclear absorption and final state interaction with 
CO- moving hadrons on the J /ij} and tp' suppression in proton- nucleus and nucleus-nucleus 
collisions. We show that a reasonable description of the experimental data can be achieved 
with theoretically meaningful values of the cross-sections involved and without introducing 
any discontinuity in the J /ij} or tj}' survival probabilities. 
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In 1986 Matsui and Satz [1] proposed J/V' suppression in heavy ion collisions as a 
signal of quark gluon plasma (QGP) formation. This suppression results from Debye 
screening in a medium of deconfined quarks and gluons. Shortly afterwards, the NA38 
collaboration tested this idea and found that in 0-U and S-U collisions the ratio J/ip 
over di-muon continuum decreases with increasing centrality [2] . However, very soon, two 
alternative explanations involving conventional physics - i.e. no phase transition - were 
proposed. In the first one, known as nuclear absorption [3] [4], the cc pair wave packet 
produced inside the nucleus, is modified by nuclear collisions in such a way that it does not 
project into J/i/j but into open charm. It was shown that the absorptive cross-section aabs 
needed to explain the ^-dependence of J/ip production in pA collisions also does explain 
the suppression found in nucleus-nucleus. This cross-section turns out to be about 6 mb. 
In the second explanation, known as interaction with co-movers, the J/i/j produced outside 
the nucleus is surrounded by a dense system of hadrons (mainly pions) and converts into 
open charm due to interactions in the medium [5]. This interaction takes place at low 
energies - which, however, have to be large enough to overcome threshold effects. Reliable 
theoretical calculations of the corresponding cross-section af^ show that it increases very 
slowly with energy from threshold [6]. In view of that, the second interpretation has been 
progressively abandoned in favor of nuclear absorption [7]. 

Very recent data for Pb Pb collisions obtained by the NA50 collaboration show an 
anomalous J/ip suppression [8]. The ratio J/i/j over Drell-Yan (DY) is two times smaller 
than the extrapolation of the O-U and SU data based on nuclear absorption - the statistical 
significance of this discrepancy being of nine standard deviations. These data provide 
a most exciting hint of QGP formation and some interpretations in this context have 
been presented at the QM'96 conference [9]. In them a discontinuity in the J/i/j survival 
probability is assumed when some threshold of local energy density is reached. 

In this note we present an attempt to describe the observed J/ i/j suppression using the 
combined effect of nuclear absorption and interaction with co-movers - without introducing 



any discontinuity in the survival probability [10]. For the nuclear absorption we adopt the 
formalism of refs. [3, 4]. For simplicity, we use the exponential form of the J/il^ suppression 
given in ref. [4] and used in the experimental papers, namely 

Si = exp(-p L Uahs) (1) 

where p = 0.138 nucleon/fm^ is the nuclear density and L is the length of nuclear matter 
crossed by the cc pair. The more involved formalism of ref. [3] gives results very similar 
to those obtained from (1). 

We turn next to the final state interaction with produced hadrons (co-movers). A 
rigourous treatment of this interaction is very complicated and is not available in the 
literature. We use the simple treatment proposed in refs. [11]. The decrease in the spatial 
density of J/i/j at point x due to the interaction i/j-h is given by [12] 

A^ = p^{x)p^{x)at (2) 

with d'^x = rdrdyd^s, where r is the proper time, y the space-time rapidity (to be later 
on identified with the usual rapidity) and d^s an element of transverse area, af^ is the 
part of the ip-h cross-section which does not contain the J/ip in the final state, averaged 
over the momentum distribution of the colliding particles. The effect of thresholds will be 
taken into account in an effective way through the value of af^. Assuming longitudinal 
boost invariance and a dilution of the densities and p'^ of the type l/r (i.e. neglecting 
transverse expansion), we get from (2) 
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where tq is the formation time and r is the duration of the hadronic phase. 

We want to express the densities dN/dyd^s in terms of the observables dN/dy. We 
have [11] 



where b is the impact parameter of the coUision and the geometrical factor G{b) is given 

by 



G{b) = 



jd^sTl{s)Tl{b-s) 



(5) 



which has an obvious geometrical interpretation. For the nuclear profile T^ib) we use 
standard Saxon- Woods. For the proton we use a Gaussian profile with Rp = 0.6 fm. Using 
(3)- (5) we have : 
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and, for a finite time interval, 
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We have now to specify the duration time of the interaction t(6). This time is not well 
known. We use the following ansatz [11]. In the case of Gaussian profiles one has 
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where R are the rms radii. In ref. [13] it was found that the quantity in brackets in (8) is 
precisely the geometrical HBT squared transverse radius. Following the arguments of [14] 
we take it as a measure of the duration of interaction and therefore use 



r{b) = [2nG{b)]-^/'^ 



(9) 



For the formation time tq we take tq = 1 fm [14] . Of course our results depend on the value 
of Tq. However, this dependence can, to a large extent, be compensated by a small change 
of afg. Finally dN'^{b)/dy is a measurable quantity. In order to avoid model estimates we 
use the experimental value of Et as a measure of the hadronic activity. More precisely, for 
SU collisions, where the NA38 calorimeter covers a range — 1.3 < r]cm < 1-1 not far from 
the one of the dimuon (0 < 77 < 1), we take 



where Et is the average energy of neutrals measured by the calorimeter in each centrality 
bin, Ar] = 2.4, and < pt > = 0.35 GeV. Note, however, that our results depend only on the 
product afgdN'^/dy. The value of b in each bin is determined from the NA38 code [15]. 
Unfortunately, in Pb Pb collisions, the calorimeter does not have the same acceptance 
as in SU and, moreover, is not located at mid-rapidities. Due to these differences the 
Et measured in Pb Pb has to be multiplied by a factor 2.35 ± 0.15 [16] in order to be 
comparable to the Et measured in SU. Finally for a pA collision we take 

^^M^/i 3 dN^^^''~ 

where P is the average number of collisions. From (1) and (7) we obtain the combined 
result of nuclear absorption and destruction of the J/ip via interactions with co- moving 
hadrons, as 

dN'^ dN'^ 
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Note that dN'^/dy at tq is close but not identical to AB times the corresponding value in 
pp collisions. This is due to the fact that, contrary to Si, S2 ^ I for pp. Therefore it has 
to be determined from 
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We can now compute the absolute yield of J/i(^ in any reaction - or the ratio J/i(^ over 
DY since, in the latter case, = S'2 = 1. The results, which depend on two parameters 
aabs and af^, are presented in Fig. 1 and compared with the NA38 and NA50 data. 
The agreement with experiment is reasonably good. In particular the strong suppression 
between SU and Pb Pb is obtained with no discontinuity in the parameters. However, our 
L-dependence is somewhat too weak in pA collisions and too strong in SU. It is important 
that the values of the parameters 

(7^6. = 4.1 mb , at = 0.mmh (14) 



are very reasonable. Had we needed a much larger value of af^ our interpretation of the 
J/il^ suppression should be dismissed on theoretical grounds [6]. 

So far we have considered, besides nuclear absorption, all the destruction channels 
h + ip ^ D + D + X, ... , with cross-section af^. Likewise, in order to study ip' suppression 
we have to consider the channels h + ip'^D + D + X,... , which do not involve ip' in 
the final state. The corresponding cross-section will be denoted af^^ . Due to the different 
geometrical sizes, afg is larger than af^ at very high energies. Their difference is even 
bigger at low energies due to the dramatic differences in the energy behaviour of these two 
cross-sections near threshold [6] . With this sole extra parameter at our disposal, it is not 
possible to reproduce the V'VV' ratio in both SU and Pb Pb systems. If we choose af^ 
such as to reproduce the SU data, the result for central Pb Pb is an order of magnitude 
too low. However, in this case the above destruction channels are not the only relevant 
ones. One has also to consider the exchange channels 



1p + 7T >1p' + X , Ip' + TT >1p + X 



(15) 



with cross-sections af^. and af^ respectively. Asymptotically, af^ = af^. However, at 
low energies af^. is expected to be much larger than af^. due to the different thresholds. 
The presence of these channels has little effect on the J/ip over DY ratio but it changes 
considerably the ip' /ip one and allows to cure the problem mentioned above. Indeed, for 
central Pb Pb collisions, when the ip'/ip ratio becomes very small, channels (15) produce 
a feeding of ip' at the expense of ip, thereby increasing the ratio ip'/ip. 

Let us now discuss the combined effect of all destruction and exchange channels. The 
destruction channel for the ip' is treated in the same way as for the ip - with af^ replaced 
by (^to i^co > ^co)- For the exchange channels (15), there is a gain of tp' due io ip ^ ip' 
conversion and a loss of ip due to the inverse reaction. The net gain of ip' is 
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where R{b) is the ratio of ip' over ip rapidity densities at time tq. The net gain of ip is 
obviously given by the same eq. (16) with opposite sign. 
Combining (1), (6) and (16) we have 
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Contrarily to (6), eqs. (17) (18), have to be solved numerically, because, not only 
changes with increasing r, but also R{b). Therefore, it is not possible to get a close formula 
at freeze-out time r - but only the variation during an infinitesimal interval At. One has 
to solve the problem numerically, dividing the total £nT interval into a very large number 
of subintervals, and using as initial condition in each subinterval the result obtained at the 
end of the previous one. 

The results for the ratios J/ip over DY and ip'/ip are given in Figs. 1 and 2. We have 
used the following values of the parameters 

aabs = 4.1 mb , at = 0.40 mb , a^o = 2.6 mb , af^ = 0.1 mb , af^ = 0.65 mb (19) 

The result presented in Fig. 1 for the J/ip over DY ratio is not changed by the introduction 
of the exchange channels (15) (within 1 %). More precisely, a small change in af^ from 
0.46 (14) to 0.40 mb (19) has compensated for their effect. The value of af^ is basically 
determined from the data on ip'/'ip for SU. Finally, the value of af^ is determined in such 
a way to get enough feeding of ip' from ip in Pb Pb. Due to the smallness of R{b), our 
results are rather unsensitive to the ratio crf^/crf^ and, in order to decrease the number of 
parameters, we have taken it equal to crfo/^co — 6-5. (A ratio crfx/^ex — 1 with af^. — 0.06 
also gives acceptable results). Although we have not attempted a best fit of the data we 
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describe the V'VV' ratio reasonably well. In particular, we have a mild decrease of this 
ratio both in pA and Pb Pb collisions and a faster decrease in SU. This striking feature is 
also present in the experimental data. However, our V'VV' ratio in pA collisions decreases 
somewhat faster than the experimental one. 

Before concluding it should be noted that the values of Et measured in Pb Pb colli- 
sions, relative to those measured in SU, are 20 to 30 % larger than expected from scaling 
in the number of participant nucleons and from Monte Carlo codes. At present this point 
is not well understood either theoretically or experimentally. If the Et values in Pb Pb 
were to be decreased by such an amount, the values of the ratio i/j' /i/j in Pb Pb would 
increase without spoiling the agreement with experiment. However, the ratio J/i/j over 
DY for Pb Pb collisions would increase (by as much as 20 % in the most central bin of Pb 
Pb) as shown in Fig. 1. In this case, the mechanism described above would not reproduce 
entirely the NA50 data. 

Note also that an important part of the effect of the co-movers comes from the region 
of T near tq where the densities are very high and one can wonder whether such a dense 
system can be regarded as a hadronic one. In any case our mechanism of J/i/j suppression 
is different from Debye screening. 

In conclusion, combining nuclear absorption and final state interaction with co-moving 
hadrons, we have obtained a reasonable description of the J/ip and ip' data. This des- 
cription is better for SU and Pb Pb collisions than for pA. It has been achieved with 
theoretically meaningful values of the cross-sections involved and without introducing any 
discontinuity in the J/i(^ or if;' survival probabilities. 
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Figure Captions 

Fig. 1 The ratio Bn^a{J /il))/a{DY) versus the interaction length L in the final state for 
pp, pA^ SU and Ph Pb collisions. The data are from ref. [8] . The theoretical values 
are obtained from eq. (17) with the values of the parameters in (19). The same 
result (within 1 %) is obtained from eq. (12) with the values of the parameters 
in (14). The straight line corresponds to nuclear absorption alone (eq. (1)), with 
aabs = 6.2 mb. 

Fig. 2 The ratio B ,j_^a{tl)') / B ^^a{ J / ip) versus L in pp, pA, SU and Pb Pb collisions. The 
data are from ref. [8]. The theoretical values are obtained from eq. (18), with the 
values of the parameters in (19). 
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